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_~ Abstract 

The recent discovery of a SM-likc Higgs boson at the LHC, with a mass around 125-126 GeV, 
i— | together with the absence of results in the direct searches for supersymmetry, is pushing the 

rG SUSY scale (tosuSy) into the multi-TeV range. This discouraging situation from a low-energy 

SUSY point of view has its counterpart in indirect SUSY observables which present a non- 
decoupling behavior with msusY- This is the case of the one-loop lepton flavor violating Higgs 
decay rates induced by SUSY, which may remain constant or even increase as TOgusY grows, 
depending on the class of intergenerational mixing in the slepton sector which are taken into 
account (LL, LR, RL or RR). In this work we focus on the LFV decays of the three neutral 
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MSSM Higgs bosons h, H, A —> r/i, considering the four types of slepton mixing (£f 3 L , <5 Lu 
^23 L j ^23^); an d show that all the three channels could be measurable at the LHC, being h — > r/i 
the most promising one, with hundreds of events expected with the current LHC center-of- 

^ mass energy and luminosity. The most promising predictions for the future LHC stage are also 

£f) included. 
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1 Introduction 

The absence of any experimental signal, so far, in the searches for supersymmetry (SUSY) at the 
LHC [1] and the discovery of a new Higgs-like particle by ATLAS [2] and CMS [3] with a mass 
rriH SM ~ 125-126 GeV, are pushing the SUSY mass parameters above the 1-TeV range. On the one 
hand, the present lower mass bounds for squarks of the first and second generations and for gluinos 
are already above 1 TeV, and on the other hand, if the observed Higgs boson is identified with 
the lightest Higgs boson h of the Minimal Supersymmetric Standard Model (MSSM), a radiatively 
corrected mass m^ ~ 125-126 GeV also implies rather heavy squark masses of the third generation 
(mainly stop masses) at or larger than 1 TeV. In principle, to place the SUSY masses at the 
multi-TeV range seems discouraging, both from an experimental point of view due to the inability 
to detect SUSY directly, and from a theoretical point of view, in regard to the naturalness of 
the theory, which contrarily suggests a soft SUSY-breaking scale, msiJSY, at or below the TeV 
scale. However, leaving the naturalness issue aside, the MSSM scenarios with very heavy SUSY 
masses can have other interesting aspects [4]. In particular, this is the case of specific Higgs boson 
observables, like certain Higgs partial decay widths, which present a non-decoupling behavior with 
ttt-susY) as shown, for instance, in [5-9], opening a new window to the indirect detection of heavy 
SUSY. As it is well known, the decoupling of SUSY radiative corrections in the asymptotic large 
SUSY mass limit is valid for SUSY theories with an exact gauge symmetry, in agreement with 
the general decoupling behavior of heavy states in Quantum Field Theory as stated in the famous 
Appelquist-Carazzone theorem [10]. Nevertheless, it is also known that this theorem does not apply 
to theories with spontaneously broken gauge symmetries, nor with chiral fermions, which is the 
case of the MSSM. Furthermore, in order to have decoupling, the dimensionless couplings should 
not grow with the heavy masses. Otherwise, the mass suppression induced by the heavy-particle 
propagators can be compensated by the mass enhancement provided by the interaction vertices, 
with an overall non-decoupling effect, which is exactly what happens in some MSSM Higgs boson 
decays to fermions. For instance, it was studied in [5] how non decoupling appears for large SUSY 
masses in the h — > bb decay through one-loop SUSY-QCD corrections, when the involved SUSY 
particle masses grow simultaneously with a generic soft SUSY-breaking scale msusY (see also [6] ) . 
A similar non-decoupling behavior was obtained for flavor changing neutral Higgs boson decays into 
quarks of the second and third generations through both SUSY-QCD corrections [7] and SUSY-EW 
corrections [8]. Other interesting non-decoupling SUSY-EW effects have also been seen in lepton 
flavor violating (LFV) Higgs boson decays within the context of the MSSM-seesaw model [9]; and 
in Higgs-mediated LFV processes like: r — > 3/i decays [11], some semileptonic r decays [12,13] and 
in [i — e conversion in heavy nuclei [14]; all of them within the MSSM-seesaw model too. Other 
non-decoupling effects from heavy SUSY particles have also been noticed within the context of the 
SUSY inverse-seesaw model [15]. Some studies of SUSY non decoupling within the MSSM have 
also been performed in the effective field theory approach where the effective Higgs-fermion-fermion 
vertices are induced from one- loop corrections of heavy SUSY particles [6, 16]. 

In the present work, motivated by the recent discovery of the SM-like Higgs boson, we will focus 
on the study of the LFV Higgs boson decays and we will work within the context of the MSSM at 
one-loop order with the hypothesis of general slepton flavor mixing. We will extend the study to 
the three neutral Higgs bosons of the MSSM h, H and A, considering the new Higgs-like particle 
to be the lightest Higgs boson h. In particular, we will study the LFV Higgs decays h — > Tfi, 
H — )• Tfi and A — > t^jl. This kind of processes provide an unique window into new physics due to 
the strong suppression of flavor violation in the Standard Model (SM), where the flavor mixing is 
induced exclusively by the Yukawa couplings of the corresponding fermion sector. This is specially 
interesting in the lepton sector in which the flavor mixing will be hugely suppressed because of the 



smallness of the lepton Yukawa couplings. Therefore, the discovery of any process which violates 
the lepton flavor number would be a clear signal of physics beyond the SM. LFV is nowadays a 
very active field which is being studied in different models, through several channels (for a review, 
see for instance [17]): radiative decays r — > /jl-j, t — >■ e^f and \i — > e^/; 3-body decays r — > 3/j,, 
t — > 3e and \i — > 3e; /j, — e conversion in nuclei; semileptonic r decays; among others. The specific 
case of LFV Higgs decays within supersymmetric models has also deserved special attention in the 
literature [9,18]. 

Our purpose here is to take advantage of the mentioned non-decoupling behavior with msuSY 
in the LFV Higgs decay widths into charged leptons of different generations, T{<f> — > Ulj), with i ^ j 
and 4> = h,H, A, in order to look for sizeable branching ratios which can give rise to detectable 
signals at the LHC. Here and from now on, Ulj with i ^ j in the final state of the LFV decays 
refers to either Irfj or l^lj. In general, the radiatively corrected LFV Higgs couplings to leptons 
are proportional to the heaviest lepton mass involved, being this the reason why we select the 
TfjL channel as the most promising one. In addition, the lie channel is extremely restricted by 
the associated radiative decay, fi — > ej [19], leaving us almost no room to move in the allowed 
parameter space of slepton flavor mixing, and driving us to extremely low rates, not measurable 
in any LHC detector. The re channel, on the other hand, gives us very similar results to the r[i 
channel, and from an experimental point of view, the LHC sensitivity to the former should be 
equivalent to the latter [20]. Therefore the results obtained along this work for the LFV Higgs 
decays into r/i are straightforwardly translated into the re channels. For the present study we 
will focus then on the (f> —> t/j, decays within the MSSM at one loop with general slepton flavor 
mixing between the second and the third generations, without assuming any particular source 
of flavor mixing. Consequently, this general slepton mixing will be parametrized in a model- 
independent way by four dimensionless parameters, S^ , $23 , $23 and <5|| . It is important to 
remark that our calculations will be made in the mass eigenstate basis for all MSSM particles, 
including both charged sleptons and sneutrinos with the corresponding exact diagonalization of 
the full charged slepton and sneutrino mass matrices, and therefore we will not use here the Mass 
Insertion Approximation (MIA). The most important constraints to these four slepton mixings 
come from the r — >• /U7 searches [21]. A recent study on updated constraints to all these general 
slepton flavor mixings 5f, B in [22] indicates that, indeed, the present bounds to BR(r — > /ry) lead 
to constraints on the 23 mixings which for 500 GeV < w-susy < 1500 GeV and 5 < tan/3 < 60 
are at \S^\ max , \S^\ max , \S^ 3 L \ max ~ O(10" 2 - 10" 1 ), and [^[max ~ 0(10^ - 1), and for heavier 
SUSY masses lead to weaker bounds. We will check that by raising the SUSY mass scale into the 
multi-TeV range will turn into the needed relaxation of these bounds, since these LFV radiative 
decays present a decoupling behavior with the SUSY scale. This will allow us to explore a high 
SUSY mass scale region, where we will find very promising values for the LFV Higgs rates. 

Our paper is structured as follows: in Section 2 we describe the MSSM frame which we will 
use in the rest of the paper, introducing the flavor mixing and the scenarios we will work with. 
Section 3 is devoted to the LFV Higgs decay rates in the heavy SUSY particle mass limit, showing 
the non-decoupling behavior of these rates with msuSY- We dedicate Section 4 to the relevant 
numerical results of the LFV Higgs cross sections at the LHC, considering the effects induced by 
each type of delta, <5^ r. This will be done taking into account the conditions of the present and 
future phases of the LHC, in order to estimate the number of events one could expect for each LFV 
Higgs channel. We will finally close the paper with the conclusions, contained in Section 5. 



2 The MSSM with general slepton mixing 

In this paper we explore SUSY scenarios with general flavor mixing in the slepton sector, which 
have the same particle content as the MSSM. Within these scenarios, the most general hypothesis 
of flavor mixing assumes a non-diagonal mass matrix in flavor space, for both charged slepton and 
sneutrino sectors. The charged slepton mass matrix is a 6 x 6 matrix due to the six electroweak 
interaction eigenstates, Il,r with I = e,n, r. The inclusion, within the MSSM, of only three 
sneutrino eigenstates, vl with v = i/ e , v„, v T , reduces the sneutrino mass matrix to a 3 x 3 matrix. 
In the case in which slepton and sneutrino mass matrices were diagonal, we would still have a tiny 
flavor violation induced by the PMNS matrix of the neutrino sector and transmitted by the tiny 
neutrino Yukawa couplings, but we neglect it here. 

The non-diagonal entries in the 6x6 slepton matrix can be described in a model-independent 
way in terms of a set of dimensionless parameters 5^ B {A, B = L,R; i, j = 1, 2, 3), where L, R stand 
for the possible chiralities of the lepton partners and i,j indexes run over the three generations. 
These scenarios with general sfermion flavor mixing lead generally to larger LFV rates than in the 
so-called Minimal Flavor Violation (MFV) scenarios, where the mixing is induced exclusively by 
the Yukawa coupling of the corresponding fermion sector. This statement is true for both squarks 
and sleptons but it is obviously of special interest in the slepton case due to the extremely small 
size of the lepton Yukawa couplings. Hence, in the present case of slepton mixing, the 5£ s will 
provide the unique origin of LFV processes with potentially measurable rates. 

One usually starts with the non-diagonal 6x6 slepton squared mass matrix referred to the 
electroweak interaction basis, which we order here as (ej,, £ll, ?l, (1r, (Ir, tr), and write this 
matrix in terms of left- and right-handed blocks M~ (A, B = L, R), which are non-diagonal 3x3 
matrices, 
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with flavor indexes i,j = 1,2,3, corresponding to the first, second and third generations, respec- 
tively. 9\v is the weak angle, Mz is the Z gauge boson mass, (m^^m^^m^) = (m e ,m^,m T ) are 
the lepton masses, tan/3 = 1)2/ v\ with v\ = ("%?) and V2 = (^2)' the two vacuum expectation 
values of the corresponding neutral Higgs boson in the Higgs SU(2) doublets, Hi = (T-i^l %i) and 
H2 = (7-1.2 ^2)' ano - A* ^ s *^ e usual higgsino mass term. It should be noted that the non diagonality 
in flavor comes exclusively from the soft SUSY-breaking parameters, which could be non vanishing 
for i / j, namely: the masses m^ • ■ for the slepton SU(2) doublets {vu lu), the masses m^ ■■ for 

the slepton SU(2) singlets (Iri), and the trilinear couplings, A\y 

In the sneutrino sector there is, correspondingly, a one-block 3x3 mass matrix, referred to the 
(v e L,D^L,D T L) electroweak interaction basis, 

Ml = ( M? LL ) , (3) 



where 



M? LLij = m 2 Lij + Qm!cos2/?) Sij . (4) 



It should be also noted that, due to SU(2)l gauge invariance, the same soft masses m^ ■ ■ enter 
in both the slepton and sneutrino LL mass matrices. Besides, in the previous equations we have 
neglected the neutrino mass terms, which due to their extremely small value are totally irrelevant 
for the present computation. 

The general slepton flavor mixing is then introduced via the non-diagonal terms in the soft- 
breaking slepton mass matrices and trilinear coupling matrices, which are defined here as 
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In all this work, for simplicity, we are assuming that all 5^ parameters are real. Therefore, 
hermiticity of M 2 and M. 2 implies 5^- B = fiE ■ Besides, in order to avoid extremely large off- 
diagonal matrix entries we restrict ourselves to \8f- | < 1. It is worth to have in mind for the rest 
of this work that our parametrization of the off-diagonal in flavor space entries in the above mass 
matrices is purely phenomenological and does not rely on any specific assumption about the origin 
of the MSSM soft mass parameters. In particular, it should be noted that our parametrization for 
the LR and RL squared mass entries connecting different generations (i.e. for i ^ j) assumes a 
similar generic form to the LL and RR entries. For instance, M? „ = SkFm-f m~ . This implies 

to ' ILR23 Z6 L 2 E 3 r 

that our hypothesis for the trilinear off-diagonal couplings A\, with i ^ j (as derived from Eq.(6)) 
is different than other possible definitions considered in the literature. In particular, it is clearly 
different to the assumption M~ . . ~ wMstjsy ; where v 2 = v\ + v \ and Msusy is a typical SUSY 
mass scale, as it is done for instance in [16]. 

The next step is to rotate the sleptons and sneutrinos from the electroweak interaction basis to 
the physical mass eigenstate basis, 
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with R l and R u being the respective 6x6 and 3x3 unitary rotating matrices which yield the 
diagonal mass-squared matrices as follows, 
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In the numerical part of the present study we will restrict ourselves to the case where there 
are flavor mixings exclusively between the second and third slepton generations, thus we set all 



5(j B, s to zero except for ij = 23. On the one hand, the LFV one-loop corrected Higgs couplings are 
proportional to the heaviest lepton mass involved [9] and, therefore, the Higgs decay rates into fie 
are suppressed by a factor in^/m^ with respect to the h,H,A—t Tfi, re decay rates. On the other 
hand, the related LFV radiative decay fi — > e^ has a much more restrictive upper bound [19] than 
r — > ej and r — > fv) decays [21], and the present allowed values of the #12 s would not drive us to 
any measurable eft — > fie rates. 

In order to simplify our analysis, and to reduce further the number of independent parameters, 
we will focus in the following numerical study on simplified SUSY scenarios, where the relevant 
soft-mass parameters are related to a single SUSY mass scale, msusY- In particular we choose the 
following setting for the relevant mass parameters: 

m L = m E = m SUSY , (11) 

fi = M 2 = amsusv , (12) 

where a is constant coefficient that we will fix in the next sections to different values, namely, a = |, 
|, 1. We also set an approximate GUT relation for the gaugino masses: 

M 2 = 2Mi = M3/4 (13) 

Here and in the following we use a short notation for the common soft masses, namely, m^ for 
m^ = m^ = m^ = m^ , etc. For simplicity, we have also assumed vanishing soft-trilinear 
couplings in the slepton sector, i.e., A T = A^ = A e = 0. We have checked that other choices with 
non-vanishing values for any of these couplings do not alter the conclusions of this work. On the 
other hand, the values of the soft masses of the squark sector are irrelevant for LFV processes, 
except in the present case of LFV MSSM Higgs bosons decays where these parameters enter in the 
prediction of the radiatively corrected Higgs boson masses. Since we want to identify the discovered 
boson with the lightest MSSM Higgs boson, we will set these parameters to values which give a 
prediction of rrih compatible with the LHC data. Specifically, we fix them to the particular values 
tjIq = rrifj = m^ = At = Aj, = 5 TeV, which we have checked provide a value for m^ that lies 
within the LHC-favored range [121 GeV, 127 GeV] for all the MSSM parameter space considered 
here. Consequently, the input parameters of our simplified SUSY scenarios are: the mass of the 
pseudoscalar Higgs boson, tua, the ratio of the two Higgs vacuum expectation values, tan/3, the 
generic SUSY mass scale, msusY) and the four delta parameters, S^s , ^23^1 ^23 an d ^23^1 which we 
vary within the following intervals: 

• 200 GeV < m A < 1000 GeV, 

• 1 < tan/3 < 60, 

• 0.5 TeV < msusv < 10 TeV, 

• "I < ^, *§?, * ^ < I- 

3 Results for the LFV decay rates 

In this section we study the behavior of the radiative corrections from SUSY loops to the LFV 
neutral Higgs bosons decays in the heavy SUSY particle mass limit and compare them with the 
case of LFV radiative lepton decays. For the forthcoming estimates of the LFV Higgs decay rates, 
we use the complete one-loop formulas and the full set of diagrams contributing to the T(<j> —> Ulj) 
and T{(j) — > lih) partial decay widths, with i ^ j, within the MSSM, which are written in terms 



of the mass eigenvalues for all the involved MSSM sparticles, including the physical slepton and 
sneutrino masses, mr (i = 1, ..,6) and m^ (i = 1,2,3), and the rotation matrices R l and R u of 
Eqs. (9) and (10). We take these general one-loop formulas from [9] and emphasize that they are 
valid for the general slepton mixing case considered here, with all the mixings effects from the S^- B, s 
being transmitted to the LFV Higgs decay rates via the physical slepton and sneutrino masses and 
their corresponding rotations. The off-diagonal trilinear couplings in Eq.(6) also enter into this 
computation of the LFV Higgs decay rates. For a more detailed discussion about these analytical 
results we refer the reader to [9]. It should be also noted that, since we are assuming real <5^ s, s, 
the predictions for BR(<^> — >• lilj) and the CP-conjugate BR(</> — > lilj) are the same. Thus, we will 
perform our estimates for just one of them and will denote this rate generically by BR(0 — > lilj). 
Obviously, in the case that these two channels cannot be differentiated experimentally one should 
then add the two contributions to the total final number of events. However, for the present 
computation we do not sum them and report instead the results for <j> — > lilj, meaning that they 
are valid for any of the two cases. 

As said above, we focus only on h — > Tfi, H — )• t/j, and A — >■ t[l decay channels and consider 
the constraints imposed over the parameter space by the current upper bound on the related LFV 
radiative decay r — > ^7 [21]. The SUSY mass spectra are computed numerically with the code 
SPheno [23]. The slepton and sneutrino spectra are computed from an additional subroutine that 
we have implemented into SPheno in order to include our parametrization of slepton mixing given 
by the ^^'s. The LFV decay rates are computed with our private FORTRAN code in which we 
have implemented the complete one-loop formulas for the LFV partial Higgs decay widths of [9] and 
the complete one-loop formulas for the LFV radiative r decay widths which we take from [11]. Note 
that these latter formulas for the r — >■ /try decays are also written in terms of the physical sparticle 
eigenvalues and eigenstates and do not neglect any of the lepton masses. The mass spectrum of 
the MSSM Higgs sector, with two-loop corrections included, and their corresponding total widths 
are calculated by means of the code FeynHiggs [24]. 

Next we present the numerical results for the branching ratios of the LFV decays. We show 
in Figure 1 the behavior of the branching ratios for the two types of LFV decays, BR(</> — > t/j,) 
and BR(r — > ^7), as functions of msiJSY; and we consider two different values of tan/3, namely, 
tan/3 = 5 (left panels) and tan/3 = 40 (right panels). In each case we set one single delta to be 
non vanishing with the particular values: 5^ = 0.5 (upper panels), S^ = 0.5 (middle panels) and 
^23 = 0-5 (lower panels). All the other flavor changing deltas are set to zero. We find identical 
results for 8^ = 0.5 as for 8^ = 0.5 and, for brevity, we have omitted the plots for 8^ in Figure 1. 

On the upper panels of Figure 1 , when the responsible for the flavor mixing between the second 
and the third generations is 8^ , the branching ratios of the LFV Higgs decays show a clear non- 
decoupling behavior with m-suSY) which remain constant from msuSY — 1 TeV to tosusy = 10 
TeV, with values of BR(/i ->r|i)~5x lO" 11 and BR(H, A->r/i)~8x 10~ 9 for tan/3 = 5, and 
BR(/i -)■ TfJ,) ~ 3 x 10~ 9 and BR(H, A ->• t/j,) ~ 3 x 10~ 6 for tan /3 = 40. Another important feature 
that should be noted is the fast growing with tan /3 of these decays which increase the H and A 
LFV decay rates almost three orders of magnitude from tan /3 = 5 to tan /3 = 40. Furthermore, we 
have numerically checked that for large values of tan /3 > 10 the partial decay widths T(H, A — > t/j,) 
go approximately as ~ (tan/3) [9], whereas T(h — > Tfi) goes as ~ (tan/3) 2 . This implies that the 
corresponding branching ratios go all at large tan/3 > 10 as BR(/i, H, A — > Tfi) ~ (tan/3) 2 in this 
LL case, since the total widths go as T tot (H, A) ~ (tan/3) 2 and r to t(^) is approximately constant 
with tan/3. In contrast, the branching ratio of the r —> fij decay presents a decoupling behavior 
with msuSY, decreasing as ~ V m susY' an< ^ ^ ^ s reduced around five orders of magnitude from 
^"-SUSY = 500 GeV to msusY = 10 TeV. In all these figures we have also included, for comparison, 
the experimental upper bound for the r — > fij channel, whose present value is BR(r — >• ^7) 
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Figure 1: BR(/i — > Tfj,), BR(H — > Tfi), BR(^4 — > r/x) and BR(r — > ^7) as functions of msusY for 
tan/3 = 5 (left panels) and tan/3 = 40 (right panels) with S^ = 0.5 (upper panels), 5%$ = 0-5 
(middle panels) and 5^3 = 0.5 (lower panels). The results for 5^ = 0.5 (not shown) are identical 
to those of 5^ = 0.5. In each case, the other flavor changing deltas are set to zero. In all panels, 
rriA = 800 GeV and the other MSSM parameters are set to the values reported in the text, with 
M2 = gW-susY- The horizontal dashed line denotes the current experimental upper bound for 
r -» ^7 channel, BR(r -)• ^7) < 4.4 x 10~ 8 [21]. 



< 4.4 x 10 8 [21]. Thus, each tosusy point which leads to a prediction of BR(r — >• ^7) above this 



line is excluded by data. Therefore, only values of msusY <] 2 TeV for tan /3 = 5 and msusY > 5 TeV 
for tan/3 = 40, and their corresponding predictions for the LFV rates, are allowed for 5^ = 0.5. 

On the middle panels of Figure 1 we have plotted the LFV Higgs and radiative decay rates as 
functions of m-susY; considering 5^ as the responsible for r — \x mixing. A similar non-decoupling 
behavior to the #§3 case can be observed for 5^ ; , whose branching ratios stay again constant 
as mgusY grows. However, the numerical contribution of S^ to the LFV processes is much less 
important than that of S^, and all the RR rates are in comparison around two orders of magnitude 
smaller than the LL rates. In the RR case, another important feature is that all the predictions 
found of BR(r — > /17) for 5 < tan/3 < 40 and msusY values above 500 GeV are allowed by the 
present BR(r — > /17) experimental upper bound. 

The predictions of BR(/i — > 771), BR(H — > r/x), BR(j4 — > Tfi) and BR(r — > /X7) as functions 
of msusY) for the case 5^ = 0.5, are shown on the lower panels of Figure 1. We see clearly 
that the contribution of #23 to the LFV Higgs decays is more relevant than the b^ and 5§^ 
ones. Firstly, for the same numerical input mixing value of 0.5, one obtains much larger LFV 
Higgs decay rates for the LR case, in several orders of magnitude, than for the LL and RR cases. 
Secondly, the most relevant feature of these LFV Higgs rates with LR mixing is their manifest 
non-decoupling behavior with msusY, which leads to an increase of the decay rates as msusY 
grows, in contrast to the LL and RR cases with constant decay rates. On the other hand, the 
BR(t — > /X7) rates display again a decoupling behavior, but the decreasing with msusY in this LR 
case is slower than in the previous LL and RR cases, going now as ~ l/ m !usY- We nave a l so 
studied numerically the behavior with tan/3 of the LFV rates in this LR case. We have found 
that for large values of tan/3 > 10, the partial widths go as F(H,A — > Tfi) ~ (tan/3) 2 whereas 
T(h — > Tfi) is approximately constant with tan/3. This implies that the corresponding branching 
ratios BR(/i, H, A — > r/x) stay all approximately constant with tan/3 at large tan/3 > 10. We have 
also checked that the radiative r decay rates are also constant with tan /3 [22] . Thus, from now on, 
we will take advantage of these two complementary effects (non decoupling of LFV Higgs decays 
and decoupling of LFV radiative decays) in order to obtain sizeable LFV Higgs branching ratios 
whilst the r —> fij rates stay below its experimental upper bound. For example, on the lower left 
panel of Figure 1, corresponding to tan/3 = 5, BR(r — > fi^) is below its experimental upper bound 
for msusY > 7 TeV and BR(/i — > Tfi), BR(H — > r/u) and BR(j4 — > r/x) can reach large values, 
allowed by data, up to 5 x 10~ 5 , 1.5 x 10~ 5 and 6 x 10~ 5 , respectively. If we increase to 40 the 
value of tan/3 (lower right panel), BR(r — > /X7) falls down its upper bound for rragusY ^ 8 TeV. 
Thus, for tan/3 = 40 and msusY = 10 TeV, we obtain allowed values of BR(/i — > r/x) ~ 4 x 10~ 5 , 
BR(H -> r/x) ~ 1.2 x 10" 4 and BR(A ^ r/x) ~ 5 x 10" 4 . 

To sum up the main results of Figure 1, the most relevant (5^ parameter is 5^ (and £23 )> 
which gives rise to much larger LFV Higgs decay rates than 5^ and ^23^. In addition, if one wants 
to obtain sizeable and allowed by data branching ratios, one needs large values of mgusY> which 
plays a double role: on the one hand, it increases the values of the LFV Higgs decay rates (due to 
the non-decoupling behavior of these decays with mgusv) an d, on the other hand, it brings down 
r — > /X7 below its experimental upper bound (because of the decoupling effect of LFV radiative 
decays with msuSY)- 

Now, we are interested in investigating if other choices of M2 alter these previous results. For 
this purpose, we have plotted in Figure 2 the predictions of BR(/i — > r/x) (dots in upper panels), 
BR(H — > Tfi) (dots in medium panels), BR(^4 — > r/i) (dots in lower panels) and BR(r — > /X7) 
(crosses in all panels) as functions of msusY for different values of a (see Eq. (12)), a = 1 (in red), 
a = 3 (in green) and 0=5 (in blue), with 5^ = 0.5 (left panels) and 5^ = 0.5 (right panels). In 
the case of LL mixing, all the LFV Higgs rates, which present the same behavior with a, increase 
around a factor of 7 from a = h to a = 1, while r — > ^7 rates present the opposite behavior with 
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Figure 2: LFV Higgs decay rates (dots) and BR(r — > /ry) (crosses) as functions of ttt-susy with 
^23 = 0-5 (left panels) and 5^ = 0.5 (right panels) for different choices of M2: M2 = msuSY (in 
red), M2 = ^msvsY (in green) and M2 = gmgusY (i n blue). The results for b§i L = 0.5 (not shown) 
are identical to those of #23 = 0.5. In each case, the other flavor changing deltas are set to zero. 
In both panels, tra = 800 GeV, tan/3 = 40 and the other MSSM parameters are set to the values 
reported in the text. The horizontal dashed line denotes the current experimental upper bound for 
r ->■ fij channel, BR(r ->■ ^7) < 4.4 x 10~ 8 [21]. 



a, decreasing in a factor about 40 for the same values of a. Therefore, if b^k is the responsible for 





Figure 3: BR(h — > r/i), BR(# — > t/j,), BR(A -> t/j.) and BR(r -> /ry) as functions of 5^ (upper 
left panel), S§^ (upper right panel), 5^ (lower left panel) and <5|| L (lower right panel). In each 
case, the other flavor changing deltas are set to zero. In all panels, tha = 800 GeV, tan/3 = 40, 
"^SUSY = 10 TeV and the other MSSM parameters are set to the values reported in the text, 
with M2 = TTtsuSY- The horizontal dashed line denotes the current experimental upper bound for 
r ->■ /ry channel, BR(r —- /ry) < 4.4 x 10~ 8 [21]. 



the slepton mixing, and for the explored interval 1/5 < a < 1, the larger M2 is (and consequently 
Mi and //), the larger the LFV Higgs branching ratios are and the lower BR(r — > fij) is. On the 
right panels of Figure 2, the results for BR(/i, H, A — > r/i) and BR(r — > fij) in the Li?-mixing case 
are depicted. Again BR(/i, H,A^ Tfi) rise as a grows and the enhancement in this LR case is 
by a factor of about 30 by changing a = g to a = 1. Different to the LL case, BR(r — > ^7) also 
increases with a, although softer than the LFV Higgs rates and, for example, for M2 = msusY and 
^2? = 0-5, all the predicted rates for this LFV radiative decay are above its experimental upper 
bound. We learn from Figure 2 that the best choice, for a fixed delta parameter, in order to obtain 
the largest LFV Higgs rates is M2 = w-susy- However, we must be very careful, especially in the 
LR case, because it is possible that these large rates are excluded by the r — > ^7 upper bound, 
depending basically on the specific values of S^ and tan/3. 

In order to look into the largest values of 5^ allowed by data for the choice M2 = tbsusy, we 
show in Figure 3 the results of the branching fractions for the LFV Higgs decays into r/i and the 
related LFV radiative decay r — > fij as functions of the four deltas considered along this work. 
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For completeness, we have also presented here the results for 5^, which are irrelevant for the 
present work since all the branching ratios obtained are extremely small to be detectable at the 
LHC, and the 5^ results, which are identical to the 5^ ones. The plots in Figure 3 show the 
expected growing of the LFV rates with the |<523 I 's, and are all clearly symmetric 5^ — > — #23 ■ 
On the upper left panel we have the results for LL case and it is clear that all the values of 
<5^ 3 L , from — 1 to 1, are permitted by data, due to the large suppression that r — > ^7 suffers for 
"^SUSY = 10 TeV. The predictions for H — >• Tfi (green crosses) are indistinguishable from A — >■ Tfi 
ones (blue asterisks), which are superimposed in this plot. One can reach values of BR(/i — > t/j,) 
~ 1(T 7 and BR(#, A 4t/i)~2x 10" 4 at the most for 8%f = ±1. The predictions for the LFV 
rates as functions of 8^ are presented on the lower left panel of Figure 3. In this case, values 
of J #23* J > 0.3 are not allowed by the r — > ^7 upper bound and the largest value of <5|^ ~ 0.3 
allowed by data gives rise to a branching fraction slightly lower than 10 -3 for the h — > Tfi channel, 
while BR(H, i-fr/j)-5x 1(T 3 . 

To close this section, we can conclude from Figures 1, 2 and 3 that, for the explored intervals of 
the model parameters, the largest LFV Higgs rates that are allowed by the r — > ^7 upper bound 
are obtained for the following set of parameters: rrasusY = 10 TeV, M2 = msusY and 8^ = 0.3. 
In the following this will be our favorite reference point and it will be explored further in the next 
sections. The corresponding predictions for other choices of #23 > -^2> an d wisiiSY can be easily 
derived from Figures 1, 2 and 3. Similarly, in the LL case we will fix 8^ = 0.9 as our reference 
point from now on. The corresponding predictions for other choices of #23 can be easily derived 
from Figure 3. 

4 Results for the LFV cross sections and event rates at the LHC 

In this section we present the results of the LFV cross sections and event rates at the LHC which 
are mediated by the production of neutral MSSM Higgs bosons and their subsequent LFV decay 
into T[i. The production cross sections of the neutral Higgs bosons are calculated here by means 
of the code FeynHiggs [24]. For low values of tan/3, the production cross sections of the three 
neutral MSSM Higgs bosons are dominated by gluon fusion. By contrast, for moderate and large 
values of tan f3 (> 10), the production cross section of H and A Higgs bosons associated with a pair 
of bottom-antibottom quarks become the dominant ones, while the h production cross section is 
still dominated by gluon fusion. In the following, we consider center-of-mass energies at the LHC 
of y^s = 8 TeV and -y/s = 14 TeV and focus in the two cases with the largest LFV Higgs decay 
rates, with either LL or LR slepton r — /i mixing. Although we do not expect any competitive 
background to these singular LFV signals at the LHC, a more realistic and devoted study of the 
potential backgrounds should be done, but this is beyond the scope of this paper. 

4.1 Predictions for single LL mixing 

We first consider the case of single LL mixing which is the delta parameter best motivated from 
theory. In fact, most of the top-down approaches to LFV physics, where the mixing parameters are 
generated through RGE running, as for instance in SUSY-seesaw models, predict that the largest 
mixing comes from ^23 , with the following hierarchy: 8^ ^> #23 3> ^23 • This motivates us to 
investigate further this case and to look for potential sizeable branching ratios and measurable 
cross sections at the LHC with just this single 5^ \ ■ The best numbers which we have been able 
to achieve for #23 are ) however, not very encouraging, as displayed in Figure 4, which corresponds 
to the following choice of parameters: M2 = wisusy, ^23" = 0-9, tan/3 = 60 and rriA = 1000 GeV 
(in order to be in accord with CMS searches [25]). On the left panel of Figure 4 one can observe 
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Figure 4: Left panel: BR(/i — > t/j,), BK(H — > r/i), BR(A — > Tfi) and BR(r — > ^7) as functions 
of msusY- The horizontal dashed line denotes the current experimental upper bound for r — > ^7 
channel, BR(r — > ^7) < 4.4 x 10~ 8 [21]. Right panel: a{<fi) X BR(</> — >• r/x) as a function of msusY> 
where <\> = h,H,A, for y^s = 8 TeV (dots) and yfs = 14 TeV (crosses). The production cross 
section of the light Higgs boson h is calculated via gluon fusion and the H and A production cross 
sections are calculated associated with a pair of bottom-antibottom quarks. The shaded gray area 
is excluded by the r — > /ry upper bound. In both panels, itia = 1000 GeV, tan/3 = 60, <5^ 3 L = 0.9 
and the other MSSM parameters are set to the values reported in the text, with M2 = m-susY- 



again the non-decoupling behavior of LFV Higgs rates with tosusy and the opposite decoupling 
behavior of the LFV radiative decay channel, which allows us to be safe of the r — > ^ upper 
bound for values of m-susY ^ 3.8 TeV. In this case, the largest rates which we obtain are yet too 
small to be detected at the present LHC phase, BR(/i — > r/i) ~ 5 x 10~ 8 and BR(H , A — > T[i) 
~ 3 x 10~ 4 . On the right panel of Figure 4 we show the cross sections for all the h,H,A — > Tfi 
channels, for the current LHC phase (dots) and for the future one (crosses), where the shaded gray 
area is the excluded region by the r — > fj.'y upper bound. As we can clearly see in this plot, the low 
cross sections predicted in the allowed region of msusY indicate that total integrated luminosities 
of C > 333 fb _1 and C > 33 fb _1 would be required at least for the H, A — >■ 171 channels, in the 
current LHC phase and in the future one, respectively, in order to obtain some event if the unique 
source of lepton flavor violation is 623 • For the h — > t\i channel it would be needed even larger 
values of C, namely, 1 x 10 3 fb _1 in the present stage of the LHC and 400 fb _1 in the future one. 



4.2 Predictions for single LR mixing 

Next we consider the case of single LR mixing that, as we will see, is much more promising than 
the previous LL mixing case. The predictions of the cross sections of LFV Higgs processes at the 
LHC with ^fs = 8 TeV are exhibited in Figure 5, for m A = 800 GeV, M 2 = msusv, 8^ = 0.3, 
tan/3 = 5 (left panel) and tan/3 = 40 (right panel). For tan/3 = 5, cr(gg — > h) x BR(/i — > r/i) can 
reach values up to 10-15 fb for large values of msusY ~ 10 TeV, due to the enhancement of the 
LFV Higgs rates as msuSY grows. The predictions for a(gg — > H,A) x BR(H, A — > Tfi) also rise 
with mgusYi and reach maximum values of about 10~ 3 fb and 2 x 10~ 3 fb, respectively, since the 
value of mA considered here is rather heavy (800 GeV) and consequently leads to small production 
cross sections for both H and A Higgs bosons. For tan/3 = 40, the prediction of the light Higgs 
boson h does not change barely but the H and A cross sections enlarge more than one order of 
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8 TeV. Right panel: a(<f> 
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boson h is calculated via gluon fusion and the H and 
associated with a pair of bottom-antibottom quarks, 
MSSM parameters are set to the values reported in the 
are excluded by the r — > /ry upper bound [21]. 



function of msuSY) where <fi = h,H,A, for 

gs production cross sections are calculated 

x BR(0 — > t/jl) as a function of msuSY, 

production cross section of the light Higgs 

A production cross sections are calculated 

for yfs = 8 TeV. In both panels the other 

text, with M.2 = msuSY- Shaded gray areas 



magnitude, with a maximum value of 7 x 10~ 2 fb for mgusY ~ 10 TeV, since both the production 
associated with bb and also BK(H,A — > r/x) increase with tan/3, as we have shown in Figure 1. In 
both panels, the shaded gray areas correspond to values of BR(r — > z/7) above its present upper 
bound and consequently excluded. 

We estimate the number of expected events at the LHC for h,H,A—> t/j, channels as functions 
of m-susY hi Figure 6. As in the previous Figure 5, the shaded gray area is excluded by the r — >• fij 
upper bound. The current phase of the LHC, with a center-of-mass energy of \fs = 8 TeV and 
a total integrated luminosity of C = 25 fb -1 (dots), and the future one, with y/s = 14 TeV and 
C = 100 fb -1 (crosses), are considered, both with the following choice of parameters: mA = 800 
GeV, tan/3 = 40, 5^ = 0.3 and M2 = tisusy- In the former phase, one is able to achieve, for the 
allowed large values of m-susY, a few hundred of expected events for the h —> Tfi channel, reaching 
~ 300 events for msuSY = 10 TeV. For the H,A^-th channels we could expect 1 or 2 events at 
the most in the current LHC phase. In the next phase of the LHC, the expectations for the number 
of events will increase hugely. Our predictions in this case (crosses in Figure 6) are between 2 x 10 3 
and 3 x 10 3 events for h — > t/jl and between 30 and 60 events for the H,A—^t/jl channels, for the 
allowed region of msuSY and the same choice of parameters. 

The low number of events obtained for the H, A — > t^jl channels are due to the large value of 
rriA considered up to now, which reduces the production cross sections of these two heavy Higgs 
bosons. Now we are going to scan over a wider region of m A values, from 200 GeV to 1000 GeV, 
fixing msusY to our reference largest value of 10 TeV and setting S^ = 0.3 and M2 = mgusY- 
Figure 7 shows the results of <t(4>) x BR(0 — > t^l) as a function of m A , where <f> = h,H,A, for 
tan/3 = 5 (left panel) and tan/3 = 40 (right panel), with a LHC center-of-mass energy of yfs = 8 
TeV. As in Figure 5, for low values of tan/3, the Higgs production cross sections are dominated by 
gluon fusion, whilst for large values of tan /3 the production cross section associated with a pair of 
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Figure 6: Number of events of (ft — > r/x as a function of m-susY; where (ft = h,H, A, for tha = 
800 GeV, tan/3 = 40 and 5^3 = 0.3. The production cross section of the light Higgs boson h is 
calculated via gluon fusion and the H and A production cross sections are calculated associated 
with a pair of bottom-antibottom quarks, for yfs = 8 TeV and C = 25 fb -1 (dots) and for y/s = 
14 TeV and C = 100 fb -1 (crosses). The SUSY parameters are chosen as in Figure 5. The shaded 
gray area is excluded by the r — > ^7 upper bound [21]. 
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Figure 7: Left panel: c(gg —> 4>) X BR(0 — >■ r/u) as a function of niA, where (ft = h,H,A, for 
"^SUSY = 10 TeV, tan/3 = 5 and 5^ = 0.3. The Higgs production cross sections are calculated via 
gluon fusion for y/s = 8 TeV. Right panel: a {(ft) x BR(c/> — > t/j,) as a function of tua, for msuSY = 
10 TeV, tan/3 = 40 and 5^ = 0.3. The production cross section of the light Higgs boson h is 
calculated via gluon fusion and the H and A production cross sections are calculated associated 
with a pair of bottom-antibottom quarks, for yfs = 8 TeV. The SUSY parameters are chosen as in 
Figure 5. The shaded blue areas are excluded by CMS searches for MSSM neutral Higgs bosons 
decaying to tT pairs in the ttia — tan /3 parameter space for the MSSM m™ ax scenario [25] . 
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Figure 8: Left panel: number of events of (j> — ¥ t/i as a function of tha, where <f> = h,H,A, for 
^SUSY = 10 TeV, tan/3 = 5 and 5^ = 0.3. The Higgs production cross sections are calculated via 
gluon fusion. Right panel: number of events of 4> — ► Tfj, as a function of tha, for tostjsy = 10 TeV, 
tan/3 = 40 and 5^ = 0.3. The production cross section of the light Higgs boson h is calculated 
via gluon fusion and the H and A production cross sections are calculated associated with a pair 
of bottom-antibottom quarks. In both panels, the number of LFV events are calculated for yfs = 
8 TeV and C = 25 fb" 1 (dots) and for ^ = 14 TeV and £ = 100 fb" 1 (crosses). The SUSY 
parameters are chosen as in Figure 5. The shaded blue areas are excluded by CMS searches [25]. 



bb quarks is the relevant one for H and A Higgs bosons. We can observe in both panels that the 
cross section of the LFV h channel is practically independent of mA, since the corresponding light 
Higgs mass stays within the interval rrih = 125 ± 4 GeV, and leads to practically constant values of 
a(gg — > h) x BR(/i — > t/j,) around 10-15 fb. However, the cross sections of LFV H and A channels 
are strongly dependent on niA, as expected, and fall down more than three orders of magnitude 
from niA = 200 GeV to tua = 1000 GeV. For low values of tan/3 (left panel), the largest cross 
section is, by far, for the h — > tjjl channel, even for low values of tha- The cusp appearing around 
m,A — 350 GeV for the A — > t^jl channel is due to the well known threshold effect when crossing by 
two on-shell top quarks. This effect also occurs for H production via gluon fusion, but it is softer, 
resulting in a plateau instead of a cusp. For the LFV A channel, we reach cross sections of 1 fb at 
the most, around the top-quark threshold effect, while the largest LFV H cross sections, for values 
of tua < 250 GeV, are excluded by CMS searches for MSSM neutral Higgs bosons decaying to tt 
pairs in the tua — tan/3 parameter space for the MSSM m™ ax scenario [25] (shaded blue area). For 
large values of tan /3 and low values of tua ;$ 300 GeV, the largest cross sections are obtained for 
the H,A—¥t/j, channels, as can be seen on the right panel of Figure 7. Nevertheless, the choice of 
parameters which give rise to these large cross sections are excluded by CMS searches [25], and for 
tan/3 = 40 only values of tua > 800 GeV are allowed by data. In this case, the cross sections of 
H and A Higgs bosons are very suppressed, as we have seen in Figure 5, and become again much 
smaller than the light Higgs boson cross section. 

Figure 8 is devoted to the number of events expected for <fi — > r/x channels as functions of tjia, 
where <ft = h,H, A. The choice of parameters is the same as in Figure 7, for tan/3 = 5 (left panel) 
or tan/3 = 40 (right panel), and we consider the present phase of the LHC (dots), with y/s = 8 
TeV and C = 25 fb -1 , and the future one (crosses) too, with \fs = 14 TeV and C = 100 fb -1 . The 
predictions of the h — >■ tjjl channel are very similar for the two considered values of tan /3, obtaining 
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around 300 events in the present LHC phase and more than 3 x 10 3 in the future one. On the 
other hand, the expectations for the H,A—>Tfi channels depend strongly on tan/3. For tan/3 = 5, 
the largest numbers for A — > T[i are achieved around the ti production threshold, reaching up to 
30 events in the present LHC phase and almost 300 events in the future phase. For the H —> Tfi 
channel, we obtain at the most 6 (70) in the present (future) phase of LHC for tua — 250 GeV, 
which is the lowest value of itia allowed by CMS searches with tan/3 = 5 [25]. These numbers for 
the heavy Higgs bosons LFV decays are hugely reduced as tua grows, due to the large suppression 
on the production cross sections. For tan/3 = 40 and low values of tua < 300 GeV, one could 
obtain more events for the heavy Higgs bosons than for the light one, due to the enhancement with 
large tan/3 in the H and A production cross sections associated with bb. However, values of tua < 
800 GeV with tan/3 = 40 are excluded by CMS searches [25] (shaded blue area in the plot). For 
the allowed tua range, the H,A—>Tfj, channels give rise to less than 2 events in the present stage 
of the LHC. If we consider the future LHC phase, this situation is enormously improved, expecting 
around 50 (20) events for the H,A-t rfi channels with itia = 800 GeV (1000 GeV). 

4.3 Predictions in the (ttia, tan/3) and (tosusy, tan/3) planes 

Once we have set up the most relevant parameters for the present study of LFV at the LHC, which 
are the two flavor mixing deltas S^ and 5^ (and correspondingly S^f 1 ), "isusy and tua masses 
and tan/3, we will present next the results for the final rates at the LHC, both in the present 
and future phases, in the most convenient way for comparison with future experimental analysis, 
namely, in the (m^,tan/3) and (msuSYjtan/3) planes. 

We start this analysis with the LL case and consider the (tua, tan/3) plane. The analysis in the 
(msusY;tan/3) plane does not provide new relevant information in this case since, as we have seen 
before, the LFV Higgs rates are nearly constant with mgusY at large msusY > 2 TeV. We have 
plotted in Figure 9 the number of events expected in the (itia, tan/3) plane for the A — > t[l channel 
with £23 = 0.9 and msusY = 5 TeV, considering both present and future LHC phases (left and 
right panels, respectively). We do not show here the results for the H — )■ rfi rates because they are 
very similar to the A — > r\i ones and do not provide any additional relevant information. On the 
other hand, the h — > t/j, channel cannot supply any significant signal at the LHC in this LL case, 
as one can infer from Figure 4, unless extremely large total integrated luminosities were considered 
(larger than 500 fb _1 ). Due to the CMS exclusion region in the (m^tan/?) plane [25] (shaded 
blue area in both panels), it is evident that we cannot expect any LFV event in the present phase 
of the LHC if the responsible for r — fi mixing is the 5^ parameter and l^ | < 1. Fortunately, 
the numbers for the future phase of the LHC are more promising, as shown on the right panel of 
Figure 9. For instance, for rriA — 450 GeV and tan/3 ~ 15, we could expect at least 1 event, and 
up to 5 for larger values of itia and values of tan/3 > 50. 

Next we analyze the results for the most promising case of LR mixing in the two (itia, tan f3) 
and (msusY) tan/?) planes. Figures 10 and 11 summarize the results for the h — >■ t\i and A — ► r\i 
channels, respectively, in the present LHC stage. Figures 12 and 13 present the results for the same 
LFV Higgs channels in the future LHC stage. Again, the corresponding results for the H — )■ t/j, 
channel are not shown here because they are very similar to the A — > t/j, channel and do not provide 
additional information. The Higgs boson production cross sections are calculated again via gluon 
fusion for h and in association with bb for A. As in the previous figures, the shaded blue areas are 
excluded by the recent CMS searches for MSSM neutral Higgs bosons decaying to tt pairs in the 
so-called m™ ax scenario [25] and the shaded gray areas are excluded by the present r —> /ry upper 
bound [21]. 

On the left panel of Figure 10, where the number of expected events from the h —> T[i channel 
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Figure 9: Number of expected LFV events in the (mA,tan/3) plane from A — > t\x for S^ = 0.9 
and m-susY = 5 TeV. Left panel: present phase of the LHC with -y/s = 8 TeV and C = 25 fb _1 . 
Right panel: future phase of the LHC with -^/s = 14 TeV and C = 100 fb _1 . In both panels the 
other MSSM parameters are set to the values reported in the text, with Mi = msusY- The shaded 
blue areas are excluded by CMS searches [25] . 



in the present phase of the LHC are shown as a function of ttia and tan/3, for msusY = 10 TeV 
and 5^ = 0.3, we see that the maximum allowed number of events are obtained in the low tan/3 
region. Hundreds of events are expected, up to 700 for tan/3 (< 6), independently of the value of 
vtla- In any case, in all the allowed region the number of predicted events are softly dependent on 
rriA and at least 270 events are obtained, even for large values of tan/3 and txia- The results in the 
(msusY;tan/3) plane, which are summarized on the right panel of Figure 10 for txia = 600 GeV 
and the same value of ^ = 0.3, show that the number of LFV expected events for the allowed 
region by data, are nearly independent on tan/3, for tan/3 > 10, and we would obtain at least 200 
and up to 400 events coming from this h — > r[i channel. 

The corresponding results for the A ^ t[x channel, displayed in Figure 11, show a very different 
behavior with tua, m-susY and tan/3 than the previous h case. The number of expected LFV events 
at the LHC via A — >• t/j, decays, N(A — > r/i), diminish as m^ increases (see left panel), due to 
the suppression in the production cross section of a heavy pseudoscalar Higgs boson; and raise 
with tan /3, due mainly to the growing of the A production cross section associated with a pair of 
bottom quarks and also of BR(^4 — > Tfi) with this parameter, as previously illustrated in Figures 7 
and 1, respectively. On the other hand, N(A — > T[i) also increases with msusY) because of the 
non-decoupling effect of LFV Higgs decays with this mass parameter. For very low values of ttia, 
around 200 GeV, we can expect up to 15 events (see left panel) in the allowed region by data (in 
white). These event rates could be hugely increased for large values of tan /3, but these are excluded 
by recent CMS searches [25]. For m^ = 600 GeV (right panel), the minimum value of tan/3 giving 
rise to some event is tan/3 ~ 12, with msusY — 10 TeV. For larger values of tan/3, from 15 to 30, 
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Figure 10: Number of expected LFV events in the present phase of the LHC from h — > Tfi for ^ = 
0.3, y^ = 8 TeV and £ = 25 fb . Left panel: results in the (tua, tan/3) plane for msusY = 10 
TeV. Right panel: results in the (msuSY> tan/3) plane for tua = 600 GeV. In both panels the other 
MSSM parameters are set to the values reported in the text, with M^ = msusY- The shaded blue 
areas are excluded by CMS searches [25] and the shaded gray areas are excluded by the present 
r — > ^7 upper bound [21]. 



one could obtain at the most between 2 and 5 events. 

Figures 12 and 13 show the same calculations as in Figures 10 and 11, respectively, but increasing 
the center-of-mass energy and luminosity to the values of the future phase of the LHC, with y/s = 
14 TeV and £ = 100 fb -1 . It should be noted that the allowed regions in the two panels remain the 
same as in Figure 10. In Figure 12 the predictions for the h — )• t^l channel show the same behavior 
with respect to the two pair of parameters as in Figure 10 but with an increase in the number of 
events of around one order of magnitude. On the left panel we obtain again the maximum amount 
of events for the lowest tan/3 values, being these nearly independent on tua-, and the rates decrease 
as we raise tan/3, showing a small variation with respect to txia for the allowed region by data (in 
white), as in the previous mentioned plot. Specifically, we obtain up to 5 x 10 3 events for tan/3 ~ 
2, and between 4 x 10 3 and 2.75 x 10 3 events for the region between tan/3 = 6 and tan/3 = 60. On 
the right panel we observe again the nearly independence with tan /3 and the increase with msusY 
as on the right panel of Figure 10. The number of events increase from 50 to 4 x 10 3 as we go from 
1 to 10 TeV in msusY, although values of tbstjsy ^ 8.5 TeV are excluded by the r — >• /ry upper 
bound. 

Similarly, Figure 13, containing the predictions for the A —> t/j, channel in the future LHC 
phase, shows an analogous behavior to that of Figure 11 in the two studied planes. On the left 
panel, the number of expected events decrease from 250 to 0.1 for low tan/3 = 4 as we go from 
m a = 200 GeV to tua = 1000 GeV. In the high tan/3 region (i.e. for tan/3 = 50) the decrease 
could go from 9000 to 25 for the same values, but since wia is excluded below 800 GeV, we obtain a 
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Figure 11: Number of expected LFV events in the present phase of the LHC from A — > t\i for 
#23 = 0.3, y/s = 8 TeV and C = 25 fb _1 . Left panel: results in the (rriA, tan {3) plane for msusY = 
10 TeV. Right panel: results in the (msijSY,tan/3) plane for rriA = 600 GeV. In both panels the 
other MSSM parameters are set to the values reported in the text, with M2 = msusY- The shaded 
blue areas are excluded by CMS searches [25] and the shaded gray areas are excluded by the present 
r — > ^7 upper bound [21]. 



maximum value of 50 events for this region. Regarding the results in the (msuSY) tan/3) plane, on 
the right panel, the non-decoupling effect is again observed, and the amount of LFV events raise 
with msusY- In this case the results are much more promising than in the present phase of the 
LHC case, obtaining a relevant number of events even in the low tan/3 allowed region. For instance, 
for tan/3 ~ 30 and msusY = 10 TeV, one could obtain 125 events. 

5 Conclusions 

The discovery of a new Higgs-like particle at the LHC is concentrating a lot of efforts in studying 
its properties, couplings and decays, in order to investigate if there is new physics behind it. In 
that sense, lepton flavor violating Higgs decays are of special interest, since they would clearly 
imply the existence of physics beyond the SM. We have discussed in this paper the possibility 
of obtaining sizeable LFV Higgs rates, induced by heavy SUSY, and detectable at the LHC. In 
particular, we have studied in detail the most interesting LFV decays of the three neutral MSSM 
Higgs bosons: h — > Tfj,, H — >• t/j, and A — > t\i. We have shown that these three channels present a 
non-decoupling behavior with msusY- On the one hand, if e)^ or ^ are the responsible for the 
intergenerational mixing in the slepton sector, these LFV Higgs decay rates remain constant with 
^SUSY- On the other hand, and even more interesting, when the slepton mixing is produced by 
%z ( or ^23^)' these decay rates do increase, instead, as msusY grows. In contrast, the related LFV 
radiative decay, r — > /X7, manifests a fast decoupling behavior with msusY and its rates are very 
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Figure 12: Number of expected LFV events in the future phase of the LHC from h —> t/j, for 5^ = 
0.3, i/s = 14 TeV and £ = 100 fb _1 . Left panel: results in the (rriA, tan/3) plane for msusY = 10 
TeV. Right panel: results in the (mgusYj tan/3) plane for tua = 600 GeV. In both panels the other 
MSSM parameters are set to the values reported in the text, with Mi = msusY- The shaded blue 
areas are excluded by CMS searches [25] and the shaded gray areas are excluded by the present 
r — > ^7 upper bound [21]. 



suppressed for large values of the SUSY scale. In this work, we have taken advantage of these two 
opposite behaviors with msusY hi order to reach large LFV Higgs branching ratios which are yet 
allowed by the present r — > ^7 upper bound. 

From our detailed analysis of the LFV Higgs decays, we have learned that the single RR flavor 
t — fj, mixing in the slepton sector cannot by itself provide sufficiently large rates which could be 
measurable at the LHC. The situation ameliorates slightly if we consider the LL mixing, whose 
best predictions are up to 5 LFV events for the H,A—>Tfj, channels at the future phase of the 
LHC. The most relevant delta parameter is clearly #23 (and #23 ) an d gives rise to very large and 
allowed by data h, H, A — >■ t/j, rates, for sufficiently large values of msusY> close to 10 TeV, where 
the LFV radiative r — > jjrj rates are suppressed below its present experimental upper bound. In 
this Li?-mixing case, we are able to obtain a large amount of LFV events at the LHC, being h — > Tfi 
the most promising channel. For this process we expect between 200 and 700 events in the present 
phase of the LHC. Regarding the H,A —¥ t/j, channels, 15 events are expected at the most, for 
very low values of itia close to 200 GeV and tan /3 close to 4. The situation improves noticeably if 
one considers the future phase of the LHC. In this case, with a total integrated luminosity of 100 
fb _1 , we predict thousands of LFV events from the lightest Higgs boson decay into t/j,, as much as 
5 x 10 3 events for very large msusY close to 10 TeV. For the LFV heavy Higgs bosons decays, the 
expectations are lower, but also increase in comparison with the present LHC phase, and up to a 
few hundred of events could be obtained depending on the values of tua, tan/3 and msusY- 

The impressive results presented along this work strongly suggest that a dedicated search for the 
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Figure 13: Number of expected LFV events in the future phase of the LHC from A — >• t/j, for 5^ = 
0.3, -y/s = 14 TeV and £ = 100 fb _1 . Left panel: results in the (rriA, tan/3) plane for msusY = 10 
TeV. Right panel: results in the (mgusYj tan/3) plane for tua = 600 GeV. In both panels the other 
MSSM parameters are set to the values reported in the text, with Mi = msusY- The shaded blue 
areas are excluded by CMS searches [25] and the shaded gray areas are excluded by the present 
r — > ^7 upper bound [21]. 



proposed LFV Higgs decays is extremely worthwhile and we believe that it should be further studied 
by the experiments at the LHC. If the SUSY mass scale is too heavy, as the present experiments 
are pointing out, and the SUSY particles cannot be directly reachable at the present or next future 
LHC energies, our proposal for LFV Higgs decays could provide an unique window to explore new 
physics and to find some hint of very heavy SUSY at the LHC. 
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